** Figure 2: Evolution of the Number of mission stations of the 3 denominations

** Presbyterian, Methodists, Catholic
* Input data: church_locations_panel.dta constructed from Blue Book data; Code and original data available at 
* Meier Zu Selhausen, Felix; Moradi, Alexander; Jedwab, Remi (2021), “Georeferenced Data of Christian Mission Stations, Ghana (1752-1932) ”, Mendeley Data, V1, doi: 10.17632/5r256b72dd.1


use "church_locations_panel.dta", replace

** Figure showing the number of mission stations 
drop if exist!=1
collapse (sum) exist, by(religion3 year)

reshape wide exist, i(year) j(religion3) string
rename existCatholic cath
label var cath "Catholics"
rename existMethodist meth
label var meth "Methodists"
rename existPresbyterian presb
label var presb "Presbyterians"

foreach var of varlist cath meth presb {
replace `var'=0 if `var'==.
}

sort year

twoway (line presb year if year>=1847, lcol(black) lpattern(solid) lwidth(thick)) ///
(line meth year if year>=1847, lcol(gs9) lpattern(solid) lwidth(thick)) ///
(line cath year if year>=1880, lcol(black) lpattern(dash) lwidth(thick)), ///
 leg(col(3) ) ///
	xlabel(1850 (10) 1932)  ylabel(0(100)700) xtitle("") ytitle("N Mission Stations", height(5))


graph export "Figure2.pdf", replace
